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ABSTRACT 

We perform a series of two-dimensional magnetohydrodynamic simulations of 
supernova cores. Since the distributions of the angular momentum and the mag- 
netic fields of strongly magnetized stars are quite uncertain, we systematically 
change the combinations of the strength of the angular momentum, the rotations 
law, the degree of differential rotation, and the profiles of the magnetic fields 
to construct the initial conditions. By so doing, we estimate how the rotation- 
induced anisotropic neutrino heating are affected by the strong magnetic fields 
through parity-violating effects and first investigate how the north-south asym- 
metry of the neutrino heating in a strongly magnetized supernova core could be. 
As for the microphysics, we employ a realistic equation of state based on the 
relativistic mean field theory and take into account electron captures and the 
neutrino transport via the neutrino leakage scheme. With these computations, 
we find that the parity-violating corrections reduce < 0.5% of the neutrino heat- 
ing rate than that without the magnetic fields in the vicinity of the north pole 
of a star, on the other hand, enhance about < 0.5% in the vicinity of the south 
pole. If the global asymmetry of the neutrino heating in the both of the poles 
develops in the later phases, the newly born neutron star might be kicked toward 
the north pole in the subsequent time. 

Subject headings: supernovae: collapse, rotation — magnetars: pulsars, magnetic 
field 
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1. Introduction 

Neutron stars arc created in the aftermath of the gravitational core collapse of massive 
stars at the end of their lives. After the supernova explosions, the nascent neutron stars 
receive the rapid rotations and large magnetic fields, which are believed to be observed as 
pulsars. Thus, rotations and magnetic fields should be naturally taken into account in order 
not only to clarify the explosion mechanism, but also to explain the observed properties of 
neutron stars. Here it is noted that the number of studies including such multidimensional 
aspects is increasing recently (Akiyama et al. 2003; Buras et al. 2003; Kotake, Yamada, 
& Sato 2003; Kotake et al. 2004; Ott et al. 2003; Kotake, Yamada, & Sato 2003; Yamada 
& Sawai 2003; Kotake et al. 2004; Fryer & Warren 2004; Miiller et al. 2004; Yamada, 
Kotake, & Yamasaki 2004; Liebendorfer et al. 2004). This may reflect the fact that only 
the neutrino heating mechanism assuming spherical symmetry seems difficult to produce the 
explosions (Rampp & Janka 2000; Liebendorfer et al. 2001; Thompson, Burrows, & Pinto 
2003; Liebendorfer et al. 2003). 

Among the unresolved mysteries in context of core-collapse supernovae, the physical 
origin of the pulsar kicks has been long controversial. Recent analyses of individual pulsar 
motions and observations of remnant associations between supernovae and pulsars indicate 
that supernovae explosions are asymmetric and that the neutron stars receive large kick 
velocities at birth (typically 300 — 400 km/s (Lyne & Lorimer 1994; Lorimer, Bailes, & 
Harrison 1997), with highest values greater than 1000 km/s (Arzoumanian, Chernoff, & 
Cordes 2002)). The existence of pulsar kicks are also supported by the evolutionary studies of 
neutron star and black hole binaries (Fryer & Kalogera 1997; Wex, Kalogera, & Kramer 2000; 
Mirabel et al. 2002) and by the detections of geodetic and orbital plane precessions in some 
binary pulsars (Cordes, Wasserman, & Blaskiewicz 1990; Kapsi et al. 1996). The recent X- 
ray observations have shown the correlation between the direction of pulsar motions and the 
spin axis of their supernovae in Vela and Crab pulsars (Helfand, Gotthelf, & Halpern 2001; 
Pavlov et al. 2001). However it is statistically uncertain whether the spin-kick ahgnment is 
a generic feature of all pulsars. 

So far two major classes of mechanism for these kicks have been suggested (see Lai, 
Chernoff, & Cordes (2001); Lai (2004) for a review). Any models or mechanisms should 
explain how the explosions become asymmetric and how the generated asymmetric explosions 
are related to the formation of the pulsar kicks. The first class of the mechanisms rehes on 
the asymmetric explosion as a result of the global convective instabilities (Burrows & Hayes 
1996; Goldreich et al. 1996) caused by the pre-collapse density inhomogeneities (Bazan & 
Arnett 1998) or the local convective instabilities (Janka & Mueller 1994) formed after the 
onset of core- collapse. More recently, Scheck et al. (2004) reported in their two dimensional 
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simulations that in the long duration explosions (more than a second after bounce), the 
neutrino-driven convections behind the expanding shock can lead the global asymmetries, 
which accelerates the remnant neutron star to a several hundreds of km/s (see, also Fryer 
(2004) for three dimensional calculations). In their simulations, the central protoneutron 
star is fixed as the inner boundary condition. Fully self-consistent calculations are required 
to confirm whether this assumption is valid or not. 

Here we pay attention to the second class of models in which the pulsar kicks arise 
from an asymmetry of neutrinos induced by the strong magnetic field {B > 10^''"^ G) in 
the collapsed supernova core. While the dipolar magnetic fields of most radio pulsars are 
observed to lie in the range of 10^^ ~ 10^^ G, several recent observations imply that the soft 
gamma-ray repeaters and anomalous X-ray pulsars have the strong magnetic fields as high as 
B ~ 10^5 G (Zhang & Harding 2000; Guseinov et al. 2003). So far, much theoretical efforts 
have been paid to formulate the neutrino opacity in such a strong magnetic fields (Horowitz 
& Li (1998); Arras & Lai (1999), Arras & Lai (1999)). On the other hand, it seems that little 
attention has been paid to investigate how the asymmetric neutrino heatings work globally 
in the supernova core based on the numerical simulations. In addition, most of the preceding 
studies relies on the asymmetry of the neutrinos only by the magnetic fields (see, however, 
Janka & Raffelt (1999); Duan & Qian (2004) for the simple models of atmosphere in the 
supernova core). It is noted that rotation should also contribute to produce the asymmetry 
of the neutrino heating (Kotake, Yamada, & Sato 2003). To our knowledge, the degree of 
the neutrino heating asymmetry under the presence of the rotation and the strong magnetic 
fields at the same time has not been investigated so far. This situation leads us to investigate 
how the rotation-induced anisotropic neutrino heating (Kotake, Yamada, & Sato 2003) are 
affected by the strong magnetic fields through parity-violating effects. In this study, we 
perform a scries of two-dimensional simulations of magnetorotational core collapse. In order 
to sec the parity-violating effects, we pay particular attention to the models, which are 
assumed to have the strong poloidal magnetic fields (~ 10^^ G) prior to core-collapse. Since 
the profiles of the angular momentum and the magnetic fields of such strongly magnetized 
stars are still quite uncertain (see, however, Heger et al. (2003)), we change the strength 
of the rotation and the magnetic fields systematically. Based on the simulations from the 
onset of the core-collapse through the bounce to the shock-stall, we demonstrate how large 
the anisotropy of neutrino heating could be globally in the supernova core. Furthermore, we 
hope to speculate their possible effects on the pulsar kicks. 

We describe the numerical methods in section 2. The main numerical results are shown 
in section 3. Discussion are given in section 4. We conclude this paper in section 5. 
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2. Numerical Methods and Initial Models 

The numerical method for magnetohydrodynamic (MHD) computations employed in 
this paper is based on the ZEUS-2D code (Stone & Norman 1992). The ZEUS-2D is an 
Eulerian code based on the finite-difference method and employs an artificial viscosity of 
von Neumann and Richtmyer to capture shocks. The time evolution of magnetic field is 
solved by the induction equation : dB/dt = V x (v x B), with v, B being velocity and the 
magnetic field, respectively. In so doing, the code utilizes the so-called constrained transport 
(CT) method, which ensures that the numerically evolved magnetic fields are divergence-free 
(V ■ B — 0) at all times. Furthermore, the method of characteristics (MOC) is implemented 
to propagate accurately all modes of MHD waves. The self-gravity is managed by solving the 
Poisson equation with the incomplete Cholesky decomposition conjugate gradient method. 
In all the computations, spherical coordinates are used and one quadrant of the meridian 
section is covered with 300 (r) x 50 (9) mesh points. We have made several major changes to 
the base code to incorporate the microphysics. First, we added an equation for the electron 
fraction to treat electron captures and neutrino transport by the so-called leakage scheme 
(Epstein & Pethick 1981; Bludman et al. 1982; Van Riper & Lattimer 1981; Van Riper 1982). 
The calculation of electron fraction is done separately from the main hydrodynamic step in 
an operator-splitting manner. Second, we implemented a tabulated equation of state based 
on the relativistic mean field theory (Shen et al. 1998) instead of the ideal gas EOS assumed 
in the original code. For a more detailed description of the methods, see Kotake, Yamada, 
& Sato (2003). 



2.1. Initial Conditions 

Recently, Heger et al. (2003) performed the most up-to-date stellar evolution calcula- 
tions including the effects of rotation, mixing, transport of the angular momentum, and the 
magnetic torques. Their models show that the toroidal magnetic fields are much stronger 
than the poloidal ones prior to core collapse. In our previous paper (Kotake et al. 2004), 
we considered such models as the initial conditions. However, if such initial configurations 
are true, no parity- violating effects of the magnetic fields on the neutrino heating can be 
obtained in the two-dimensional axisymmetric simulations. Considering the uncertainties 
in Heger's evolution calculations that they are one-dimensional models, which hinder the 
accurate treatment of the transport of the angular momentum and the generation of the 
magnetic fields, it should be still important to take the poloidal magnetic fields as the ini- 
tial conditions. Thus in this study, we choose to assume that only the poloidal magnetic 
fields exist prior to core-collapse, and then, prefer a parametric approach to construct the 
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initial conditions varying the profiles of the initial poloidal magnetic fields. We assume the 
following two rotation laws and the four configurations of the magnetic fields. 



1. shell-type rotation 



a;(r) = Wo X (1) 



and shell- type magnetic fields, 



Rq^ 



where a;(r) is an angular velocity, B{r) is the poloidal component of the magnetic fields, r 
is radius, and cuq, Rq and Bq are model constants. 



2. cylindrical rotation: 



V 2 '74 



u;{X, Z)^u;oX ^^^^ • i^^^, (3) 



and cylindrical magnetic fields, 

X2 + Xq^ + z, 



y- 1 74 



Z) = So X ■ ^W^' (4) 



and uniform cylindrical magnetic fields, 

= So, (5) 

where is the poloidal magnetic fields parallel to the rotational iZ) axis, X and Z denote 
distances from the rotational axis and the equatorial plane and Xo, Zq are model constants. 
The other parameters have the same meanings as above. 

Finally, we take into account the quadrupole magnetic fields following the prescription 
by Ardeljan, Bisnovatyi-Kogan, & Moiseenko (1998), 

B^{X,Z) = F^{0.5X,0.5Z -2.5) - E^{0.5X,0.5Z + 2.5), 

B^{X,Z) = F,(0.5X, 0.5Z- 2.5) -F^(0.5X, 0.5^ + 2.5), (6) 



with 



and 



^0,quad -^Ojquad 



, 2XZ 2X'^Z 
F,{X,Z) = Bo\ 



(Z2 + 1)3 (Z2 + 1)5 

FJX,Z) = Bo{^^ \ (8) 
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where is the poloidal magnetic fields parallel to the equatorial plane and Xo,quad, -^o.quad 
are model constants. 

We have computed 11 models changing the combination of the strength of the angular 
momentum, the rotation law, the degree of differential rotation, and the profiles of the 
magnetic fields. For all the models, we fix the values of £'in/|VF| to be 0.1%, where £'in/|VF| 
represents the ratio of the magnetic to the gravitational energy. The model parameters are 
presented in Table 1. The models are named after this combination, with the first letter 
,"S (shell)", "C (cylindrical)", denoting the rotation profiles, the second letter ,"L (long), 
S (short)", indicating the values of Rq and Xq, the third letter "(U) uniform, (NU) non 
uniform, (Q) quadrupole" indicating the configurations of the initial poloidal magnetic field, 
the fourth number -1 or -2, representing that the initial T/|l^| is 10~^, 10~^%, respectively. 
Note that the ratio of rotational energy to gravitational energy is designated as T/|VF|. 

We note that our choice of the initial value of T/|iy| is much smaller than the past 
magnetorotational studies in context of core-collapse supernovae (LeBlanc Sz Wilson 1970; 
Bisnovatyi-Kogan et al. 1976; Miiller & Hillebrandt 1979; Symbalisty 1984; Ardeljan et al. 
2000) to reconcile with the recent stellar evolution calculations of rotating stars (Heger et 
al. 2000). However it is still larger than the most slowest rotating models predicted by the 
most recent stellar calculations (Heger et al. 2003). As for the configurations of the initial 
magnetic fields, we prepared a variety of them in order to cover as wide the parametric 
space of the field configurations as possible. As for the initial strength of the magnetic 
fields, we pay particular attention to the models, which are assumed to have the strong 
poloidal magnetic fields (~ 10^^ G) prior to core-collapse in order to see the parity-violating 
effects. Astrophysically, such models are hkely to associate with the formation of the so-called 
magnetars. 



3. Results 

3.1. Magnetohydrodynamic Features 

We briefly summarize some dynamical aspects which will be helpful in the later discus- 
sions. For the more detailed magnetohydrodynamic computations during core-collapse, see 
also our accompanying paper (Takiwaki et al. 2004). 

In Figure 1, the entropy distributions at the final states for some representative models 
are presented. The final time of all the models is about 50 ms after core bounce (see Table 1). 
At the final time, the shock wave produced at core bounce tends to stall in the core except 
for the models, which have the centrally condensed magnetic flelds and rapid rotation with 



Table 1. The Model Parameters. 



Model 


Rotation profiles 


B profiles 


Wo (s 1) 


Bo (G) 


Ro, Xo, Zo xlO** (cm) 


SL(U)-2 


Shell 


Uniform Cylinder 


0.6 


1.9 X 10^2 


Ro = 1 


SS(U)-2 


Shell 


Uniform Cylinder 


9.0 


1.9 X 10^2 


Ro = 0.1 


CS(U)-2 


Cylinder 


Uniform Cylinder 


6.3 


1.9 X 10^2 


Xo = 0.1,Zo = 1 


SS(NU)-2 


Shell 


Shell 


9.0 


1.1 X 10" 


Ro = 0.1 


CS{NU)-2 


Cylinder 


Cylinder 


6.3 


4.7 X 10" 


Xo = 0.1,Zo = 1 


SL(Q)-2 


Shell 


Quadrupole 


9.0 


1.4 X 10" 


Ro = 1, Xo,quad = ^O.quad = 1-6 


CS(Q)-2 


Cylinder 


Quadrupole 


6.3 


1.4 X 10" 


Xo = 0.1, Zo = 1, Xo,quad = -^O.quad = 1-6 


SL(U)-1 


Shell 


Uniform Cylinder 


1.7 


1.9 X 10^2 


Ro = 1 


SS(U)-1 


Shell 


Uniform Cylinder 


28.3 


1.9 X 10^2 


Ro = 0.1 


SL(NU)-1 


Shell 


Shell 


1.7 


5.2 X 10^2 


Ro = 1 


SS(NU)-1 


Shell 


Shell 


28.3 


1.1 X 10" 


Ro = 0.1 



Note. — Note in the table that "B profiles" indicates the profiles of the initial magnetic fields. Note also that the number (-1 
or -2), which each model has in its name, indicates the values of initial T/|W| of 10"^, 10~^%, respectively. 



Table 2. Properties of the Final States. 



Model t{ (ms) Bp (G) B^ (G) (km) r^^ (km) T/\W\{(%) E^/\W\{{%) 



SL(U)-2 


40.2 


1.5 


X 


1016 


2.4 


X 


1015 


157 


175 


2.5 


X 10" 


-1 


1.9 


X 


10" 


-1 


SS(U)-2 


41.8 


3.4 


X 


1016 


2.3 


X 


1015 


178 


233 


2.5 


X 10- 


-1 


1.7 


X 


10- 


-1 


CS(U)-2 


40.7 


3.2 


X 


1016 


3.7 


X 


1015 


144 


612 


2.5 


X 10" 


-1 


1.7 


X 


10" 


-1 


SS(NU)-2 


42.5 


5.2 


X 


1016 


1.9 


X 


1015 


163 


148 


1.1 


X 10" 


-1 


5.3 


X 


10" 


-1 


CS(NU)-2 


43.1 


4.6 


X 


1016 


4.9 


X 


10" 


176 


1471 


1.1 


X 10" 


-1 


3.7 


X 


10" 


-1 


SL(Q)-2 


37.1 


1.2 


X 


1016 


9.4 


X 


1015 


146 


125 


2.6 


X 10" 


-1 


5.2 


X 


10" 


-2 


CS(Q)-2 


46.0 


3.0 


X 


1016 


9.3 


X 


1015 


169 


215 


3.4 


X 10- 


-1 


7.1 


X 


10- 


-2 


SL(U)-1 


42.3 


3.0 


X 


1016 


6.0 


X 


1015 


226 


573 




1.4 




3.8 


X 


10- 


-1 


SS(U)-1 


47.1 


4.1 


X 


1016 


2.0 


X 


1016 


133 


1361 




1.8 




3.6 


X 


10" 


-1 


SL(NU)-1 


45.0 


2.7 


X 


1016 


4.9 


X 


1015 


278 


444 




1.0 




8.3 


X 


10- 


-1 


SS(NU)-1 


44.4 


5.8 


X 


1016 


7.0 


X 


10" 


425 


1325 


2.8 


X 10" 


-1 


8.6 


X 


10- 


-1 



Note. — tf represents the final time measured from core bounce. Bp and are the maxiiiniiii strength 
of the poloidal and the toroidal magnetic fields, respectively, r^^ and are the distances from the center 
to the stalled shock front in the equatorial plane and along the rotational axis, respectively. 
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Fig. 1. — The final profiles of models SL(U)-2 (upper left), SS(Q)-2 (upper right), SS(U)-2 
(middle left), SS(NU)-2 (middle right), CS(NU)-2 (lower left), SS(NU)-1 (lower right). They 
show color-coded contour plots of entropy (ke) per nucleon. 
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Fig. 2. — The cumulative mass vs. specific angular momentum distribution in unit of Mq for 
models CS(U)-2 and CS(NU)-2 at the final state. Note in the figure that 'Initial' represents 
the initial configuration of the specific angular momentum for the models. It is seen that the 
angular momentum is transferred more efficiently for the model with the centrally condensed 
magnetic fields profile (model CS(NU)-2) than model CS(U)-2. 



Table 3. Analysis of the Heating Rate in the Strong Magnetic Fields. 



Model 


R^P (km) 


R^" (km) 


^mag (%) 


SL(U)-2 


66.6 


73.1 


0.13 


SS(U)-2 


71.4 


71.4 


0.27 


CS(U)-2 


57.3 


73.1 


0.27 


SS(NU)-2 


65.1 


68.2 


0.08 


CS(NU)-2 


54.3 


68.2 


0.18 


SL(Q)-2 


69.8 


68.2 


0.14 


CS(Q)-2 


66.6 


68.2 


0.09 


SL(U)-1 


57.3 


73.1 


0.37 


SS(U)-1 


55.8 


66.6 


0.42 


SL(NU)-1 


50.0 


65.1 


0.26 


SS(NU)-1 


51.4 


66.6 


0.12 



Note. — R„'' and Ri,^ are the distances from the 
center to the neutrino sphere, respectively, ilmag = 
\^Qt B^o/^vB=o\ maximum ratio of the 

neutrino heating rate contributed from the parity- 
violating effects, AQ^g_^Q to the heating rate without 

the corrections, QjJ'g_Q. 




angle (degree) R (cm) 

Fig. 3. — The ratio of the magnetic (pmag) to the matter pressure (pmatter) with the polar 
angle on the neutrino sphere (left panel) and the shapes of the neutrino spheres for some 
representative models (right panel). Note in the right panel that R, Z represents the distance 
from the rotational axis and the equatorial plane, respectively. 
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the strong differential rotation initially (see Figure 1 and Table 1 for models CS(NU)-2, 
SS(U)-1, and SS(NU)-l). Even if the shock waves of such models manage to approach the 
outer boundary of the numerical simulation {Z ~ 1500 km) along the rotational axis, the 
explosion energy is about one order of magnitude less than the canonical one. Thus, the 
neutrino heatings in the later phases should be necessary for the successful explosions. Here 
it is noted that the final value of T/|H^|'s and ii^m/|W^|'s for our models are well below the 
past MHD studies (Symbalisty 1984) assuming rapid rotation prior to core-collapse. Even 
in our slower rotating models than the previous simulations, it is found that the asymmetric 
matter motions can appear as seen from Figure 1. 

During core-collapse, the strength of the magnetic fields exceeds the QED critical value 
-Sqed = 4.4 X 10^^ (G), above which the neutrino reactions are affected by the parity- 
violating corrections of weak interaction. Just before core bounce, the poloidal magnetic 
fields dominate over the toroidal ones due to our assumption of the initial field configuration. 
After core bounce, the strong differential rotations near the surface of the protoneutron star 
wraps the poloidal magnetic fields to convert into the toroidal ones, whose strength becomes 
as high as -B ~ 10^^ G at the final time of the simulations. The growth of the toroidal 
magnetic fields can be estimated by an order-of-magnitude estimation as follows, 

~ 2.0 X 10- G (^) ( ^) ( ^^^^) . (9) 

where is the typical strength of the magnetic field parallel to the rotational axis just be- 
fore core bounce. At is the duration of the wrapping time, which is roughly equivalent to the 
time from core bounce to the final time, cu is the mean angular velocity at the surface of the 
protoneutron star. This estimation is roughly consistent with the numerical simulations. It 
is noted that the tangled toroidal magnetic fields {B^ < 10^^ G) are weaker than the poloidal 
ones {Bp > IQ-*^^ G) at the final time in our models (see Table 1). Magnetorotational insta- 
bility (MRI) is an another candidate to amplify the growth of the toroidal magnetic fields. 
Although the typical time scale (Balbus & Hawley 1998) corresponding to the maximum 
growth of MRI (tmri ~ 0(10) ms) at the surface of the protoneutron star is shorter than 
the simulation time scale (~ 50 ms), it is much longer than the growth time scale of the field 
wrapping (r„rap ~ 0(1) ms). Thus the amplification of the toroidal magnetic fields in our 
models is mainly as a result of the winding up original vertical poloidal fields. 

In Figure 2, the cumulative distribution of specific angular momentum is shown for 
models CS(U)-2 and CS(NU)-2. In the magnetorotational core-collapse, the anisotropic 
magnetic stress transfers the specific angular momentum of each fluid element outwards even 
in axisymmetry if the fluid rotates differentially. It is noted that the mass of the unshocked 
core is about ~ O.SMq in these models. From the flgure, it is found that the transport of the 
angular momentum occurs more efficiently for the model CS(NU)-2, whose initial magnetic 
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fields are more centrally condensed, than the model CS(U)-2 (compare the specific angular 
momentum in Figure 2). In addition, there are regions just outside the unshocked core where 
the specific angular momentum decreases outwards. In such regions, the so-called Rayleigh 
instability condition holds, which induces the convection at the later phases near the surface 
of the unshocked core in the equatorial plane. These differences mentioned above stem only 
from the degree of the initial concentration of the magnetic fields between the pair models 
and are common to the other pair of the models. 

After the final time of our simulations, we calculate the shapes of the neutrino spheres 
by the prescription (Kotakc, Yamada, & Sato 2003). Although the magnetic fields become 
strong as i? ~ 10^^ G in the vicinity of the neutrino sphere (r = 50 ~ 70 km, density ~ 
10^-^ g cm~^), the pressure of the magnetic fields (pmag = ~ 10^^(5/2 x 10^^ G)^ dyn cm~^) 
are smaller than that of the matter (pmatter ~ 10^*^ dyn cm~^) in the corresponding region 
(see the left panel of Figure 3). Thus it is found that the magnetic fields do not affect 
the shape of the neutrino sphere significantly. In fact, the neutrino spheres are deformed 
to be oblate for the rapidly rotating models, on the other hand, to be rather spherical for 
the slowly rotating models (see the right panel of Figure 3 and compare the shape of the 
neutrino sphere of model SS(NU)-1 (:faster rotating model in our computation) with that of 
the other models (:slower rotating models)). 



3.2. Effects of the strong magnetic fields on the neutrino heating 

Based on the neutrino spheres obtained in the previous section, we will estimate the 
heating rates including the corrections from the parity-violating effects outside the neutrino 
sphere. This is admittedly a very crude estimate, since it is pointed out (Liebendorfer et al. 
2001; Thompson, Burrows, & Pinto 2003) that the net heating, the heating minus cooling, 
becomes positive ~ 50 — 100 msec after the shock stagnation. In this respect, we confirmed 
that the cooling dominates over the heating also for our models including the corrections 
of the parity-violating effects. Since our calculations ended before the heating dominates 
over the cooling, it is impossible to estimate the location of the gain radius, beyond where 
the heating dominates over the cooling. Here it is noted that the contrast of the heating 
to the cooling is found to be similar in our short simulation runs. Thus, the gain radius 
should be also deformed in the later phase. If this feature lasts in the later phase, the bare 
heating rate discussed in the following and the resultant neutrino asymmetry due to the 
parity-violating corrections, will help to understand the net neutrino asymmetry. Hence, 
bearing this caveat in mind, we study the bare heating rate with the parity- violating effects 
for the final configurations (several tens msec after the bounce) in our simulations. 
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First of all, we state how we estimate the neutrino heating rate in the strongly magne- 
tized supernova core. The analysis scheme is schematically shown in Figure 4. Neutrinos are 
assumed to be emitted isotropically from each point on the neutrino sphere and stream freely 
later on. Then the heating rate at a given point outside the neutrino sphere can be found 
by summing up all the rays from the neutrino sphere. The heating by electron neutrinos 
proceeds via the charged current interaction: 

+ n — > p + e~. (10) 

We employ the corresponding cross section in a strongly magnetic field given by Arras & Lai 
(1999), which is expressed as follows : 

aB(^'^^)=(7o(^^^)( 1 + eabsn-B), (11) 

where the corresponding cross section without the magnetic fields, n, B are the 

unit momentum vector for the incoming neutrinos and the unit vector along the magnetic 
field direction (see Figure 4), and eabs is the asymmetry parameter given by 

Cabs = Cabs (e) + eabs (np) . (12) 



Here eabs(e) and eabs(np) are expressed in the leading order at the neutrino energies as 



1 hceB 



eabsiej = - 



-v 



+ 3ci' 



(13) 



and 



eabs(np) = 2 



Ca(ca + Cv) fiBnB T 



Cy + 3C^ 



kn T 



X 



Ca(ca + Cy) flBnB 



knT 



Ca(ca - Cy) flBpB 
knT 



Cy + 3C^ 



(14) 



which represent the effects of electron Landau levels in magnetic fields and polarizations of 
neutron and proton, respectively. Here k^ is the Boltzmann constant, fi^n = —1.913 fi-^, 
/iBp = 2.793 /iN ai'c the magnetic moments of neutrons and protons, respectively, with fi^ 
being the nuclear magneton, cy = 1, ca = 1.26 are the coupling constants. At the considered 
densities and temperatures, fermion phase space blocking and dense-medium effects can be 
safely ignored. The heating rate is given as 



(7b ^^^'^ c nn€^-—^d€^dn, 



(15) 
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where rin is the number density of neutron, (Priu/ {dci, (Kt) is related to the neutrino distri- 
bution in the phase space f{e^,^): 



d^n^ 1 



U^..^)4- (16) 



Introducing Eq.(16) to (15) yields 

^^"^Wfl ^^^^ 

Introducing Eq. (11) to (17) yields 

Q\ = I ao^^^'^hl Me.,rt) dn de. 



crir, 



J (To(^^^) eabs n • B Me,, Q) de, 

= Q+B=o + Agj^^o- (18) 
The first term in Eq. (18) can be expressed, 

+ _ 3ci + laocnn (el) f 3 . . x ..q^ 

where 

(el) = y" de, el Me,) ( ^ de, e^/.(e,)) ~\ (20) 

Here f2 is the solid angle, within which a point outside the neutrino sphere can receive the 
neutrinos. It is calculated by the geometric calculations (see Figure 4). As stated, the 
neutrino emission from each point on the neutrino sphere is assumed to be isotropic and 
take a Fermi distribution with a vanishing chemical potential. The Eq. (19) is expressed 
with the local neutrino flux, j,, which is in turn flxed from the neutrino luminosity obtained 
in the simulations: 



^(kT ? 

^v,B=G - 7- ^O^n , 2N2 -^^ 1^1^ 



(3 4+1) ^^ F3(0)' 
4 7r [rneC^y 
where we used the relations, 

{el)^FMlFm{kT,)\ (22) 
with ^5(0) and ^3(0) being the Fermi integrals and 



> = c -^1^ j de, el M^u) (^1^1^ 



— c 



^ Jde^elMQ- (23) 
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By the same procedure as above, the heating rate arising from the parity violating effects 
can be given as follows, 



3ci+l 



^hrrR ^V~^A ^ Ca(ca + Cv)F4(0) piM k T \ 

2^""^4T35 " ^ 4 + 3ci F3(0) ^^^^^^ ^^^^^^ 

Ca(ca - Cy) /iBp^ ^5(0) ., ^ ^2 

4 + 34 k^T F3(0)^''^'^^ 
Ca(ca - Cy) F4(0) \ 



(24) 



Here -F4(0) is the Fermi integral. The inner product of n • B can be obtained at each point 
by the geometry. 

Using the hydrodynamic quantities at the final time of the numerical simulations (~ 50 
msec) as a background, we demonstrate how much the rotation-induced anisotropic neutrino 
heating can be affected by the parity-violating corrections in the following. 

First of all, we take the model SS(U)-2 as an example model. The configuration of 
the poloidal magnetic fields at the final state is presented in the top left panel of Figure 5. 
It is seen from the panel that the poloidal magnetic fields are rather straight and parallel 
to the rotational axis in the regions near the rotational axis, on the other hand, bent in 
a complex manner, in the other central regions. This can be understood as follows. The 
core bounce occurs anisotropically due to rotation. This induces the convective motions. 
The magnetic fields are bent with the convective motions because the magnetic fields are 
frozen in fluid elements. This is because the ideal MHD is assumed in our simulations. The 
top right panel of Figure 5 shows the ratio of the neutrino heating rate corrected from the 
parity- violating effects, AQ^b-^q heating rate without the corrections, Q^^^^q. Note in 

the panel that the values of the color scale are expressed in percentage and that the central 
black region represents the region inside the neutrino sphere. It should be remembered that 
we took the two choices of T/|Ty| as the initial condition. In the models with the smaller 
T/|l^| = 10~^ %, the shapes of the neutrino sphere are found to be rather spherical and the 
resultant neutrino heating from the neutrino sphere becomes almost isotropic (see the models 
whose names end with "-2" in the right panel of Figure 3 and Table 3). In such models, it 
is found that the initial rotation rate and magnetic fleld arc not large enough to make the 
shapes of the neutrino sphere deviated from being spherically symmetric as in case of the pure 
rotation (Kotake, Yamada, & Sato 2003). In the top right panel, it is seen that the values 
of the ratio become negative in almost all the regions (see also the bottom right panel of 
Figure 5 for clarity) . This means that the heating rate is reduced by the magnetic fields than 
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Fig. 4. — Schematic drawing of the analysis scheme of the heating rate. Neutrinos are 
assumed to be emitted isotropically from each point on the neutrino sphere. The heating 
rate at P can be found by summing up the contributions of all rays, reaching P. 
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that without. Especially, this tendency is most remarkable in the regions near the rotational 
axis and the surface of the neutrino sphere (sec the regions colored by blue in the top right 
panel of Figure 5). These features can be understood as follows. Let us remember that the 
first term in Eq (24) represents the electron contribution arising from the ground landau 
level, and the other represents the contribution from the neutron and proton polarizations. 
By an order-of-magnitude estimation, it can be readily found that the first term dominates 
over the other in the central core, with typical values oi B > 10^^ G ,T ~ O(MeV). Since 
the coefficient of the first term, which is proportional to the coupling constants of weak 
interactions (oc Cy — c\, where cy — 1,ca — 1.26) is always negative, then the suppression 
or the enhancement of the heating rate through parity-violating effects is determined by the 
signs of the inner product of n • B. As mentioned above, the magnetic fields are almost 
aligned and parallel to the rotational axis in the regions near the rotational axis. With the 
rather spherical neutrino emission, the values of the inner product becomes positive in the 
regions, which results in the reduction of the heating rate. On the other hand, the values of 
the ratio is found to become positive in some regions, where the poloidal magnetic fields are 
bent due to convections (see the regions colored by red and colored by white in the top and 
bottom right panels of Figure 5, respectively). In fact, the contributions of the inner product 
are negative in such regions. In the bottom left panel of Figure 5, the contour of the ratio in 
the 360 latitudinal degrees region of a star, is prepared in order to see the global asymmetry 
of the heating. Since our simulations assume the equatorial symmetry, the above features in 
the northern part of the star become reverse for the southern part. These features on the 
global asymmetry of the neutrino heating rates are almost common to the other models. In 
the following, we state some important features of the global asymmetry produced by the 
difference of the initial rotation rate or the initial configurations of the magnetic fields. 

The configurations of the poloidal magnetic fields and the ratio of ^Q^^^q/Q^^^q for 
some representative models at the final states are shown in Figure 6. In the models, whose 
initial profile of the rotation law and the magnetic field is cylindrical with high central 
concentration (the models whose name has "CS(NU) "), the shapes of the shock wave becomes 
very coUimated near the rotational axis (see the bottom left panel of Figure 1). The final 
profile of the magnetic fields is also centrally concentrated with the strong magnetic fields 
(S ~ 4 X 10^^ G), which is favorable for the heating through parity- violating effects (see the 
top left and right panels of Figure 6). As a result, the regions, where the equatorial symmetry 
of the neutrino heating are broken, become rather wider at the north and the south poles 
than the other models, with the different initial profiles of the rotation and magnetic fields, 
while fixing the initial strength of rotation. If the initial strength of rotation is more large, 
the neutrino sphere is deformed to be more oblate due to the centrifugal force (see the black 
region of the middle right panel of Figure 6) . Since the radii of the neutrino spheres tend to 
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be smaller near the rotational axis than that in the equatorial plane (see Table 3) , the solid 
angle within which a point outside the neutrino sphere receives the rays from the neutrino 
sphere is enlarged for the vicinity of the rotational axis. This effect makes the neutrino 
heating near the rotational axis more efficient as in case of pure rotation (Kotake, Yamada, 
& Sato 2003). This rotation-induced anisotropic neutrino heating enhances the effect of 
the symmetry of the heating induced by the parity- violating corrections (see the values of 
-Rmag in Table 3 and compare the values between the models with the same name except for 
the final number, for example, SL(U)-2 and SL(U)-l). Furthermore, the regions, which are 
heated asymmetrically, are enlarged for the faster rotating model than the slower rotating 
model (compare the top right and the middle right panel of Figure 6) . 

Finally, we refer to the model whose initial profile of the magnetic field is quadrupole 
(see the bottom panels of Figure 6). In such models, it is found that there appear the regions 
with the strong magnetic fields, which seem like the horns, near the rotational axis ranging 
from ~ 100 km < Z < 200 km (see the bottom left panel of Figure 6). As a result, the 
heating regions influenced by the parity-violating corrections, are found to become rather 
off-axis (see the bottom right panel of Figure 6). 



4. Discussion 

In our computations, only electron neutrinos are treated approximately. Since all the 
numerical simulations are limited to rather early phases (~ 50 msec), it may not be a bad 
approximation. However, this is a certain limitation to the study of the later phases, which 
are more important to investigate the outcome of the north-south neutrino heating asymme- 
try discussed here. As already mentioned, if we take into the cooling, no heating region is 
found at these early times. As shown, for example, by Thompson, Burrows, & Pinto (2003), 
the heating region or (the gain region) appears more than 50 ~ 100 msec after core bounce. 
Hence, our estimate of anisotropic neutrino heating is admittedly a bold extrapolation to the 
later phases, which we have not computed yet. As stated, the multidimensional MHD sim- 
ulations which solve the neutrino transport with the parity- violating reactions are required 
in order to compute the relation between the asymmetric neutrino heating induced by the 
parity-violating corrections and the pulsar-kick. Although wc hope that the general features 
obtained in this study will be taken over later on, further numerical investigations should be 
required in order to check the validity of our results. 

As for the initial configurations of the poloidal magnetic fields, we assume that they have 
the same profile as the rotation. The assumption may not be unnatural because the angular 
momentum transfer in quasistatic stellar evolution is tightly connected to the formation of the 
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magnetic fields (Parker 1979; Spruit 2002). However, multidimensional and self-consistent 
calculations are required to be done for the understanding of the magnetorotational evolution 
of supernova progenitors. 

5. Conclusion 

We performed a series of magnetohydrodynamic simulations of supernova cores. As 
for the microphysics, we employ a realistic equation of state based on the relativistic mean 
field theory and take into account electron captures and the neutrino transport via the 
neutrino leakage scheme. Since the profiles of the angular momentum and the magnetic 
fields of strongly magnetized stars are still quite uncertain, we constructed 11 initial models 
changing the combination of the strength of the angular momentum, the rotations law, the 
degree of differential rotation, and the profiles of the magnetic fields in order to cover as 
wide a parametric freedom as possible. Based on the numerical simulations, we estimated 
how the rotation-induced anisotropic neutrino heating is changed by the parity-violating 
corrections of neutrinos under the strong magnetic fields. By so doing, we investigated the 
global asymmetry of the neutrino heating in a strongly magnetized supernova core. As a 
result, we found the following. 

1. At the prompt-shock timescale (~ 50 ms after core bounce), the toroidal magnetic 
fields are formed mainly by the wrapping of the initial poloidal magnetic fields. The strength 
of the toroidal magnetic fields is smaller than the poloidal ones at the timescale. Although 
the growth timescale of MRl corresponding to the fastest growing mode of the instability 
is within the prompt-shock timescale, it is much longer than the growth of the wrapping 
timescale. The inefficiency of the growth of MRI is due to the initial slow rotation rates 
prior to core-collapse, which reconciles with the results of the state-of-the-art stellar evolution 
calculations. 

2. Whether the parity- violating corrections in the strong magnetic fields suppress or 
enhance the heating rate than that without the corrections is found to be simply determined 
by the inner product of the unit vector of incoming neutrinos (n) and the unit vector along 
the magnetic field direction (B). In the models except those with the quadrupole magnetic 
fields initially, the poloidal magnetic fields near the rotational axis are almost aligned and 
parallel to the axis at the prompt-shock timescale. Although the neutrino emission is almost 
spherical in the slowly rotating models with the initial values of T/lVFl = 10~^ %, the values 
of the inner product becomes positive in the regions by the geometry of the magnetic field 
lines. This results in the reduction of the neutrino heating rates. In the faster rotating models 
with the initial values of T/|iy| = 10~^ %, the rotation-induced anisotropic neutrino heating 
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is found to more enhance the asymmetric heating induced by the parity- violating correction 
than in case of the slower rotating models. This is because the rotation-induced neutrino 
heating heats the regions near the rotational axis more efficiently, where the magnetic fields 
are typically strongest. The ratio of the correction to the heating rate to the heating rate 
without the magnetic fields is found to range from 0.08% to 0.42%. 

From the above results, the parity- violating corrections reduces about < 0.5% of the 
heating rate than that without the magnetic fields in the vicinity of the north pole, on the 
other hand, enhances about < 0.5% in the vicinity of the south pole of a star. Based on 
the results, we speculate that the neutron star might be kicked toward the north pole in 
the subsequent time, if the asymmetry of the neutrino heating in the north and the south 
poles develop in the later phases. This is because Scheck et al. (2004). recently pointed out 
the random velocity perturbation with an amplitude of ~ 0.1 % added artificially at several 
milliseconds after bounce can be a seed of the global asymmetry of the supernova shock 
and a cause of the pulsar-kick. In their models, the orientation of pulsar-kicks depends 
stochastically on the initial random perturbation. This is the case for the observed pulsars 
with the canonical magnetic fields {B ~ 10^^ G). However, the asymmetry of the neutrino 
heating in the strong magnetic fields {B > 10^^ G) considered here might have infiuence on 
the alignment of spin and the space velocity of the pulsar associated with the formation of 
the magnetars. Although no pulsar kicks with the very strong magnetic fields, have been 
observed so far, the alignment of spin and the space velocity of the magnetars is expected 
to be observed by future observations. 

As stated earlier, our simulations are crude in the treatment of microphysics compared 
with recent spherical models. This might be justified as long as we study the earlier phase 
in the prompt-explosion timescale. We are currently developing a two-dimensional neutrino 
transport code, which will be necessary to test the validity of implications obtained in this 
study (Kotake et al. 2004 in preparation). 
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Fig. 5. — Various quantities for the model SSU-2. Top left panel shows contour of the 
logarithm of the poloidal magnetic fields (log[-Bp (G)]) with the magnetic field lines. Top right 
panel represents the ratio of the neutrino heating rate corrected from the parity-violating 
effects, AQ^g_^Q to the heating rate without the corrections, Q^b=o- Note in the panel that 
the values of the color scale are expressed in percentage and that the central black region 
represents the region inside the neutrino sphere. Bottom right panel is the same with the 
top right panel except that the bottom right panel only shows the regions with the negative 
values of the ratio. Thus the white region shows the regions with the positive values of the 
ratio. Bottom left panel shows the contour of the ratio in the whole region, which is prepared 
in order to see the global asymmetry of the heating induced by the strong magnetic fields. 
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Fig. 6. — The contour of the logarithm of the poloidal magnetic fields (G) (left panels) and 
the contour of ^Q^b^q/Q^b=o percentage (right panels). The top, middle, and bottom 
panels, correspond to the models CS(NU)-2, SS(NU)-1, and CS(Q)-2, respectively. 



